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We analyze the properties of a Raman quantum light-atom interface in long atomic ensemble and its 
applications as a quantum memory or two-mode squeezed state generator. We consider the weak-coupling 
regime and include both Stokes and anti-Stokes scattering and the effects of Doppler broadening in buffer 
gas assuming frequent velocity-averaging collisions. We find the Green functions describing multimode 
transformation from input to output fields of photons and atomic excitations. Proper mode basis is 
found via singular value decomposition for short interaction times. It reveals that triples of modes are 
coupled by a transformation equivalent to a combination of two beamsplitters and a two-mode squeezing 
operation. We analyze the possible transformations on an example of warm rubidium-87 vapor. The model 
we present bridges the gap between the Stokes only and anti-Stokes only interactions providing simple, 
universal description in a temporally and longitudinally multimode situation. Our results also provide an 
easy way to find an evolution of the states in a Schrodinger picture thus facilitating understanding and 
design. 

Keywords: Quantum optics; Raman scattering; Quantum memory; Multimode interaction; Buffer gas; 
Warm atomic vapors; Light-atom interface 


1. Introduction 


The off-resonant Raman interface is a vividly developing approach to quantum memory An 

off-resonant anti-Stokes scattering in optically thick atomic vapor can be used to transfer excitations 
from the optical field to atoms in the read-in process j^. Atomic coherence that takes on the form 
of spinwave is created and can be converted back to l^ht in the read-out process i, i . The 
spinwaves can be stored and then further manipulated (^. A similar approach incorporates cold 
atomic ensembles instead of warm vapors . The Raman process is perhaps the simplest out of 
many possible realizations of quantum memories. In addition, a Stokes transition can be driven by 
properly tuning the pump beam, which leads to the spontaneous creation of excitations pairs in the 
write-in process, where for each photon a sing le atom is excited in a spinwave [^. The interface also 
lends itself to a multispatial mode use 10l4l2l|. Simultaneous driving of both Stokes and anti-Stokes 
transitions leads to a number of interesting phenomena enabling quantum engineering When 

they are in perfect balance, quantum non-demolition (QND) measurement type of interaction occurs 

Em. 


Here we focus on the effects of coexistence of both Stokes and anti-Stokes scattering in the 
presence of Doppler broadening. We analyze the situation for a wide range of detunings resulting 
in rather arbitrary proportions between elementary processes. It turns out that in the off-resonant, 
high-optical-depth case when losses can be neglected the combined action of the Stokes and anti- 
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Stokes scattering creates similar spatiotemporal mode structure as each of them separately would 
11,0. However, the transformation of the annihilation/creation operators of respective modes is 


more involved. It can be decomposed into a pair of beamsplitters and a squeezing operation between 
pairs of modes: atomic, photonic Stokes and photonic anti-Stokes. In particular we show that both 
read-out from atomic memory to light and read-in of light state into the memory typically involves 
significant squeezing of the quantum state. The squeezing can be only avoided when anti-Stokes 
coupling is much higher than the Stokes. The model we present facilitates understanding and design 
under realistic experimental conditions. 

This paper is organized as follows. In Section [2] we write out equations of evolution of coupled 
quantum fields. We include effects of frequent collisions with buffer gas. The solution is given 
in a form of the input-output relations for which mode basis diagonalizing the relations exists. 
We analyze the quantum evolution of the coupled triples of modes. In Section [3] we rewrite the 
interaction of three fundamental modes in a form easiest to comprehend and analyze. Section U] 
concludes the paper. 


Evolution of fields 


2.1. Maxwell-Bloch equations 

We consider a spatially extended ensemble of length L of three or more level atoms moving with 
different velocities but otherwise identical. The atoms are assumed to have two long-lived ground 
levels \g) and \h) used to store quantum information. The Raman interaction is driven by intense 
classical field P in a square-shaped pulse of duration T. We consider only forward scattering and 
thus assume no dependence of the fields on perpendicular coordinates. We also neglect all kinds of 
decoherence which is justified for large optical depths. Initially all atoms are assumed to occupy 
\g) level. As illustrated in Fig. [U^a) the interaction couples the atoms to weak quantum fields 
(sidebands) r{z,t) and corresponding to anti-Stokes and Stokes modes, respectively. 

The beams are assumed to be virtually collinear and therefore only the velocity v of the atoms 
along the beams affects the interaction. The splitting of the ground levels is taken to be so small 
that the two-photon resonance condition is simultaneously met for all velocity classes. However 
the detuning of the pump P field from resonance is different for each velocity class and therefore 
the evolution of the atomic state is somewhat different. Thus we take the atomic coherence field 
operator b^, t, v), understood within the Holstein-Primakoff approximation (0 . to be velocity class 
dependenil^. 

Together with warm atomic ensembles an inert buffer gas, such as neon, krypton or xenon is 
used to make the motion of atoms diffusive and consequently prolong the lifetime of stored spatial 


mode Hi, [2311 ■ The presence of such buffer gas inherently leads to multiple effects affecting quantum 
memories |24l |. one of which are the velocity changing collisions. These collisions transfer atoms from 
velocity class v' into v with probability per unit ti me g iven by 'jyK^v •(— v')dv', where 7 ^ is the 
collision rate and K{v v') is the collisional kernel 25l |. 

Under the assumption that the driving field P is off-resonant with given detuning from line 
centroid A, the excited level can be adiabatically eliminated. The interaction between the atoms 
and the sidebands is given by a set of Maxwell-Bloch equation s 12611 . They can be cast in terms of 
field operators in a reference frame co-moving with weak light j20||: 


df{z,t) 

dz 


j V9{v)cRiA 


+ kv)b{z, t, u)du 


( 1 ) 


^The atomic field operator is normalized so that the commutator [b{z, t, u), {z\ u')] = S{t — t')S{v — v'), therefore b{z, t, v) = 

-^==== \gi){hi\, where the summation runs through N atoms in slice Az and velocity class Au. 
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(a) 



(b) 



Detuning A [2n GHz) 

Figure 1. (color online) (a) Off-resonant Raman scattering in a 3-level atom. Linearly polarized pump of amplitude P couples 
two separate hyperfine components of the ground state, \g) and \h) via the excited state |e}. Pump is detuned from the |^} — |e} 
transition by A -|- Age and from the \h) — |e} transition by A -|- A^e? where A = 0 means the pump is tuned to the line 
centroid. Tilting and shading of the upper level symbolizes the Doppler shift and atomic population distribution respectively, 
(b) Imaginary parts of coherent scattering coefficients c/^(A), ch^(A) for rubidium 87 D1 line (A ~ 795 nm) as a function of 
detuning A between the driving field and the centroid of the line. We take two hyperfine components of the ground state as 
levels \g) = |F = 1, tuf = 0) and \h) = \F = 2, mp = 0). We take the atom number density equal N = 10^^ cm~^ and pump 
field Rabi frequency equal to the natural linewidth r/27r = 5.75 MHz. Real parts of the coherent scattering coefficients are 
found to be much smaller than the plotted imaginary parts and consequently we exclude them from consideration. Shaded 
regions close to resonances (closer than three times the Doppler width) where losses and pump absorption have significant 
effect are excluded from the analysis. 


J \/ g{v)cwA + kv)b{z, t, (2) 

^ V 9{v)cwA + kv)w^ {z, t) 

+ y/g{v)cR{A + kv)f{z, t) - s(A + kv)b(z, t, v) (3^ 

+7^, J K{v ^ v')^^^b{z,t,v')dv', 

where g{v) = ^ exp(—^■) is the thermal velocity distribution and A: = ^ is the light 

wavevector. The interaction is parametrized by coherent scattering coefficients cr{A) and cw{A) 
plotted in Fig. (Hb) and given by: 


dw\z^ t) 
dz 
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cr(A) 


cw{^) 



Pdggdgh 

(r - 2 i{A + Ahe))h 


Pdfigdeg 

(r + 2i(A + Agf^))h 


( 4 ) 

( 5 ) 


where ui = Ajj is the difference between resonant frequency of |i) — |j) transition and line 
central frequency (also see Fig. Hia) for reference), dij is the dipole moment for this transition 
and physical constants are denoted as usual. The sum runs through all possible excited states 
|e). Limiting cases correspond to \cw\ ^ |c_r| for anti-Stokes scattering and |cij| <C \cw\ for Stokes 
scattering. If the pump is far-detuned, so that A-|-Ajj ^ F, we find the imaginary parts of coherent 
scattering coefficients cr and cw to be strongly dominant over the real parts, under the assumption 
that pump amplitude P is real. Indeed, the classical field amplitude P can be always made real 
by suitable shift of the time reference. The same assumptions hold for the differential Stokes shift 
s(A) given by: 


s 


(A) = y__ 

^ ^ y 2(r - 2i(A + Ahe))h? 


+ 


\Pd 


eg I 


2(r + 2f(A +Age))^2' 


( 6 ) 


2.2. Fast collisions approximation 

The inclusion of velocity changing collision kernel in Eq. ([3]) may severely influence the solution, as 
different Maxwellian velocity distribution preserving kernels might be used (^ . 27, |2^. However, 
if we assume that the collisions are fast compared to the Raman interaction jy Lc^ ^, then 

the velocity dependence of number operator of atoms in \h) state represented by {z,t,v)b{z,t,v) 
remains close to thermal equilibrium. Consequently, we may separate out the known, Gaussian 
velocity dependence and assume b{z,t,v) = b{z, t)^yg{v). 

The equation Q for b{z,t,v) can now be integrated formally and averaged over velocity distri¬ 
bution b{z, t) = f y/g{v)b{z, t, v)dv, yielding: 


b{z, t) = b{z, 0) 

+ J (dRf{z,t')+ clYW^{z,t'fj dt', 

where the differential Stokes shift s(A) and coherent scattering coefficients cr{A) and cw(A) 
are now averaged according to the formula x = f x(A -|- kv)g{v)dv. Substituting this result to 
Eqs. ([I|)-(l2|) gives the following equation of evolution for the photonic modes: 
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+ 



e 


A ( riz,t) \ ^ 
dz \w^z,t)J 

CrC^ 

\cwcr \cw\‘^ 



6(z,0)e""‘ 


\w^z,t')) 


dt’. 


( 8 ) 


Note the average atomic coefficients cw, cr and s are virtually purely imaginary for far-detuned 
pump field. From now on we will assume these coefficients are imaginary. 


2.3. Input-output relations 

The solution of equations dZ])-® takes on a form of a linear transformation between the input and 
output quantum fields. First, a squeezing transformation i.e. hyperbolic rotation R{C) between the 
photonic modes is used to diagonalize the matrix from Eq. (|8]). Modes c\z, t) and d{z, t) are defined 
as linear combinations of f{z,t) and w^{z,t) according to the formulas: 


( r{L,t) \ 

\w^L,t)) 


R-\C) 


( c\L,t)\ 

\diL,t)J 


(9) 




f{0,t) 

t) 


where the diagonalizing hyperbolic rotation matrix is defined as 


( 10 ) 


^(C) 



sinh(C) 

cosh(C)\ 

with C = atanh 


.cosh(C) 

sinh(C )) 

for \cw\ > |cr| 


cosh(C) 

sinh(C)\ 

with C = atanh ^ 


sinh(C) 

cosh(C )j 

for Cw < c_R 


( 11 ) 


Note a significant difference between the two cases from above equation. When \cw\ > |c_r| 
both c(z,t) and d{z,t) are annihilation operators. On the contrary when \cw\ < |cj?| both c(z,t) 
and d{z,t) play role of bosonic creation operators. We can generally write [c{z,t),c^{z,t')] = 
[d{z,t),d^{z,t')] = sgn(|cw| — |cK|)(5(t — t'). In particular in the limiting cases of pure Stokes or 
pure anti-Stokes scattering corresponding to \cw\ |cr| and \c\y\ ^ | c _ r |, the R{C) is an identity 
or mode-swap transformation, respectively. 

The mode d{z,t) turns out to be decoupled and constant, thus d{L,t) = d{0,t). The equations 
coupling atoms b{z, t) with photonic mode c{z, t) are the same as in single-sideband Raman scat¬ 
tering [2 


26, 


and their solution reads: 


b{z,T)= c\0,t')y/KJ:o{z,T-t')dt' 


-|- / 6(/,0)[e ^5{z — z')-\-TT,i{z — z',T)]dz' 


( 12 ) 
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(a) write 



Figure 2. (color online) Decompositions of interaction between n-th coupled triple of modes into hyperbolic rotations R{C) 
and squeezing by exp(^n) or beamsplitter with transmission Tn- (a) Stokes domination \cw\ > \cr\ [corresponding matrix 
decomposition given by Eq. imi] and (b) anti-Stokes domination \cyy\ < |cij| [corresponding matrix decomposition given by 
Eq. lUill]. 


c^{L,t)= j b{z'— z\t)(iz' 

Jo 

+ [ c\0,t')[LJ:i{L,t-t') + 6{t-t')]dt' 

Jo 


where the interaction strength is measured by 


(13) 


K — + |cw|^ 


(14) 


and T,i{z,t) = Kzt), T,Q{z,t) = /o(2V Kzt), where In{x) are the modihed Bessel func¬ 

tions of the first kind. 


2.4. Singular modes 

The input-output relations in Eqs. m and (I13p can be simplihed by suitable choice of mode basis 
for atomic and photonic (^'[ fields in which n-th atomic mode mixes only with 

photonic modes with the same number ll|, [l^, [20|]. The annihilation operators for modes before 
and after the interaction are dehned as follows: 6^™^ = f Un’^\z)b(z,0)dz, = f Vn'^\t)r(0,t)dt, 

= f ut'^^\z)b(z,T)dz, = f vi°'^^\t)r(L,t)dt and = 

/ t)dt. The resulting operator, e.g. bn’^\ annihilates a single excitation in the re¬ 

spectful mode, in this case spatial mode Un^\z). The modes are orthonormal and the commutator 
= bmn and the same holds for all other operators. Then the relations between n-th 
input and output operators are given by a 3x3 matrix of real values equivalent to a diagram 
given in Fig. [2l 


/ bt'^ \ 

( 4 "' \ 






(15) 


If |crv| > |c_r|, the centerpiece is a squeezing operation, as in Fig. [2Ka): 
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____ n=0 n=l n=2 ____ n=3 

(a) (b) 




Position z 


Position z 


Figure 3. (color online) Atomic (a) input u'^(z) and (b) output u^^(z) (b) mode functions for first four coupled modes, 
including the fiat fundamental mode (n = 0). 


Gn = il®R{C)) 


cosh(^„) sinh(^„) 
sinh(^„) cosh(,f„) 


(-ir 


{ieR-\C)) 


(16) 


When |ci^| < |cr|, squeezing is replaced by a beamsplitter as depicted in Fig. (S^b): 


G„ = (l©i?(C)) 


\/l - Tn Tn 


-Tn 




i-iy 


{leR-Ho) 


(17) 


The mode basis for atomic {z) and photonic fields as well as squeezing or 

beamsplitter transmission for the central operation are calculated from singular value decompo- 

Normalized to length of the atomic 


sition of Green functions from Eqs. (1121) and (1131) [11 


sample L and pump pulse duration T they depend solely on the product k,LT = (c^ + |cwP)7^7"- 
For weak interaction we calculate squeezing = (kLT)”'"*'^/^ if |cvv| > |c_r| or beamsplit¬ 
ter transmission = (—if |cvk| < |ci?|- The mode functions are in the lowest or¬ 


der independent of k: 


Un^\z) = 


(in) 

Vn 


yj L(2ri+1) 


Pn{2z/L - 1), 


= Pn{l - 2z/L), 


Z = 


=p,(2 t/r-i), = 


'\J L(2n-t-l) 

v/r( 2 n+i) " ^ ~ VT( 2 n+i) ^^*'^ “ Legendre 

polynomial. We plot exemplary input and output atomic mode functions in Fig. [3l 

Note that in the limiting cases of pure Stokes or pure anti-Stokes scattering corresponding to 
\cw\ |ci?| and |cvy| <C |ci?| we recover either two-mode squeezing of bn with wh or a beamsplitter 
transformation between bn and r„. 


3. Results 

3.1. Parameters and circuit diagram 

Now we proceed with detailed interpretation of total transformation for any mode triple. Here and 
in the following we drop index n and consider only the most coupled mode n = 0. In the previous 
section we have shown the interaction of three coupled modes can be represented by squeezing or 


7 
























June 30, 2015 


Journal of Modern Optics nianuscript_jmo 


universal 



Figure 4. (color online) Universal decomposition of Raman interaction for n-th coupled triple atomic bn , photonic Stokes Wn 
and anti-Stokes Vn modes in terms of beamsplitters of transmission rjn and central squeezing by exp(xn)- Those parameters 
depend on both detuning from respective atomic resonances A and resulting ratio of the Stokes to anti-Stokes coupling cwI as 
well as the optical depth and duration of the interaction. Achievable range of parameters are plotted in Fig. [5] For simplification, 
in text parameters are given for n = 0 without index. 


beamsplitter transformation sandwiched between hyperbolic rotations as depicted in Fig. [2l The 
hyperbolic rotation can be interpreted as another squeezing therefore we effectively describe the 
interaction as a composition of squeezing operations, which makes total transformation G difficult 
to tackle. 

To facilitate understanding and design we decompose the total transformation G as a squeezing 
sandwiched in between beamsplitters depicted in Fig.|4l This form is valid for any values of cy/ and 
cr and corresponds to the following product: 




(18) 


where B{ri) 


-v 



is the beamsplitter transformation between the atomic 


mode and anti-Stokes mode with transmission ry, and S{x) 


( cosh(x) 
I sinh(x) 


sinh(x)\ 
cosh(x )) 


is the squeezing 


operation. 

The relation between parameters depicted in Fig. [2]i.e. hyperbolic rotation and squeezing ^ or 
transmission r and parameters of general transformation is found in a closed form. For the squeezing 
parameter x we obtain: 


cosh(x) 


cosh^(C) cosh(^) — sinh^(C) 
cosh^(C) — Vl — sinh^(C) 


for \cw\ > \cr 
for \cw\ < \cr 


and for the beamsplitter transmission ry: 


ry = 


l+cosh(^) 

l-l-cosh(x) 

l-hx/l — 

l-l-cosh(x) 


for \cw\ > |ck| 
for \cw\ < |ci?| 


(19) 


( 20 ) 


In Fig.[5]we plot the squeezing x and beamsplitter transmission jy as a function of pump detuning 
A which determines hyperbolic rotation ^ for various interaction strengths. We observe that not all 
transformations are achievable. An ideal readout process corresponds to using 50/50 beamsplitter 
i.e. rj = l/\/2 twice without squeezing (x = 0). At opposite extreme, spontaneous pair creation 
process corresponds to transparent beamsplitters ry = 1 around squeezing x in the first order equal 
V kLT. We see that neither is actually possible within the realistic model we present. To get as close 
as possible to pure Stokes or anti-Stokes process we need to tune the pump close to appropriate 
resonant transition, so that either \cw\ S> |cij| or |cij| 3> |cw|- 
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Detuning A (GHz) 



Figure 5. (color online) (a) Transmission of beamsplitter 77 and (b) squeezing x universal decomposition depicted in 

Fig.m and described by Eq. GHJi as a function of pump detuning A for a various interaction parameters from decomposition 
depicted in Fig. [2] squeezing ^ for domination of the Stokes \c\y\ > |c^| or transmission r for domination of the anti-Stokes 
| ch | > \cw\- 


3.2. Output states 


In this section we study the output states of quantum memory or squeezed state generator in typical 
experimental situation. First, we consider the process of spontaneous Raman scattering, ideally 
realizing two-mode squeezing. The initial state of n-th mode is |000)in regardless of mode index n. 


and \cw\ > |c_r|. The ket |A:,/,m)in/out denotes k excitations in atomic spinwave mode u 


(in/out) 




l photons in the anti-Stokes mode and m photons in the Stokes mode From 

the universal decomposition depicted in Fig. |5]we can easily obtain the output state Z//n|000)in in 
the Schrodinger picture with Un being the evolution operator. The output is most populated for 
n = 0 to which we limit ourselves dropping the index n in the following consideration. The final 
state is a two-mode squeezed exp(y) times state with one mode transmitted trough a beamsplitter 
of transmission rj: 


X 


E 


k=0 


ZY|000)in 


1 

cosh(x) 


^ tanW (x) 
j=0 





( 21 ) 


Note, that the above equation is valid for \cw\ < |cr| as well. As squeezing ^ 
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is small, we may obtain the following (unnormalized) expansion: 


Z^|000)in « |000)out +cosh(C)|101)oute 
+(cosh^(C)|202)out + cosh (C) sinh(C)|011)out)?^ + • • • 

We observe that in the first order of interaction we indeed create a pair of photon and an atomic 
excitation, which is a perfect write-in process. In the second order we not only create two such pairs, 
but also a pair of photons with each of the photons in one of the photonic modes, corresponding to 
spontaneous four-wave-mixing process with degenerate pump. Note, that the heralding of two-atom 
excitation can be made by counting the photons only in the Stokes sideband. This is possible by 
spectral filtering indeed necessary in many experiments is, El- 

Main quantum memory operation modes are read-in, that is mapping single photon in anti- 
Stokes mode to atomic mode, and the opposite read-out operation j^. Transforming the operators 
according to the known G transformation, for r <C 1 we can express the output state in the following 
way: 


^f|010)in 

2 2 

= (—rcosh(C)6^ + (1 “ ^cosh^(C))r^ + ^ cosh(C) sinh(C)tc)^^|000)in 


(23) 


Where the operators ¥, and w are bosonic creation/annihilation operators in Schrodinger 
picture acting on kets \k, I, m) in the usual way, and we calculate ZY|000)in from Eq. m- We obtain 
the following expansion for the unnormalized state: 


Z^|010)in ~ |010)out + (cosh(C)|100)out + sinh(C)|lll)out)r (24) 

As expected, no state transfer occurs in the zeroth order. In the first order photon is transferred 
to the memory, but also a pair of Stokes photon and an atomic excitation is created. Note that in 
typical experimental situations r is close to unity for many modes. In this regime contribution of 
Stokes photons might be even more significant. Very similarly, when the quantum memory contains 
a single excitation, we may calculate how the state is transformed during read-out. For small r the 
transformation is following: 


Z^|010)i 

= ((1 - rV2)it + rcosh(C)r^ — r sinh(C)rc)Z^|000)i 


(25) 


The output state can be expanded again, to show that in the first order pair of photons and 
atomic excitations are created as well: 


^I|100)i„ « |100)out + (-cosh(C)|010)out + Vcosh(2C) - l|201)out)r 


(26) 


The above calculations show that spectral filtering is extremely important, as output state will 
always be polluted by spontaneously created photons even in the lowest order for weak interaction. 
Finally, we note that the ratio of the Stokes to anti-Stokes coupling cwjcR depends solely on 


parameters of atoms. This ratio should be minimized |30|, l31[ to ensure more faithful memory 
operation, as less spontaneous Stokes photons will be emitted. 


10 




June 30, 2015 


Journal of Modern Optics nianuscript_jmo 


Conclusions 


We have introduced a model that allows us to describe Raman interaction with both Stokes and 
anti-Stokes sideband present in a long atomic sample in the presence of Doppler broadening under 
the weak interaction regime. As an exemplary medium we take warm rubidium 87 vapors contained 
in an inert buffer gas. The gas makes the motion of the atoms diffusive, allowing longer memory 
lifetime, and induces fast velocity-class thermalization, which we take into account in our model. 
Within this regime we are able to use singular value decomposition 11, [Igl . 113, HI] to find modes of 
the interaction. In this special mode bases for the field of atomic excitations, Stokes and anti-Stokes 
photons only triples of modes are coupled, one from each of the fields. We find closed expressions 
for the mode functions and coupling strengths in the weak coupling regime. 

The interaction between any triple of modes constituting Raman interface is easiest understood 
when further decomposed into a product of three two-mode quantum operations, as depicted in 
Fig. m We give closed expressions for the parameters of constituent operations. We analyzed which 
transformations are achievable in case of warm rubidium 87 vapors. 

Neither pure Stokes nor pure anti-Stokes scattering is possible, as one process is always accompa¬ 
nied by another. The effects of this interplay are easily analyzed within the framework developed. 
As an example we discuss effects on the storage in and retrieval from the quantum memory. Map¬ 
ping ideal single excitation is inevitably accompanied by creation of additional excitations due to 
four-wave mixing even when interaction is weak. These processes can be suppressed by minimizing 
the ratio of the Stokes to anti-Stokes coupling which is possible in atomic vapors to a limited extent 

We believe that our results may serve as a zeroth-order approximation for modeling various 
deleterious effects leading to decoherence such as absorption, spontaneous emission and effects due 
to motion of atoms. Multimode decoherence typically leads to quite complex results [33|. Finally 
we observe that the presented model may be applied to many similar systems, where exact atomic 
structure is different |J, 34-3 
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